Decoherence in a nonequilibrium environment 
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We consider a quantum harmonic oscillator coupled to a general nonequilibrium environment. We 
show that the decoherence factor can be expressed in terms of a measurable effective temperature, 
defined via a generalized fluctuation-dissipation relation. We further propose a simple experimental 
scheme to determine the time-dependent effective temperature in a linear Paul trap with engineered 
reservoirs. Our formalism allows quantitative description of nonequilibrium decoherence in the 
presence of an arbitrary number of non-Markovian noise sources in a unified manner. 
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Recent theoretical and experimental work on decoher- 
ence has focused on the paradigmatic model of a sin- 
gle quantum harmonic oscillator (system) linearly cou- 
pled to an ensemble of harmonic oscillators (equilib- 
rium reservoir) The interaction with the reser- 
voir leads to the dynamical suppression of interference 
fringes which arise in coherent superpositions of oscilla- 
tor states. Decoherence, therefore, plays a fundamental 
role in the description of the quantum-to-classical tran- 
sition 0, Consider for example a typical cat state, 
'^cat(^) = [^-{x) + ?/^+(x)]/V^, consisting of a linear su- 
perposition of two Gaussian wave packets {x) of width 
a and separated by a distance a). Then, for tem- 
peratures much larger than the oscillator frequency, the 
interference- fringe contrast, or decoherence factor, de- 
cays exponentially with time according to jj, 0] , 



(1) 



The coherence time tc = 1/{D(P) is here inversely pro- 
portional to the amplitude D of the thermal fluctuations 
and to the square of the spatial separation of the wave 
packets. The Gaussian dependence of the decoherence 
factor on the distance d is responsible for the nonobser- 
vation of quantum superposition of macroscropic states. 
For an equilibrium reservoir, the coefficient D is re- 
lated to the friction constant rj via the Einstein relation, 
D = TjT^ where T is the temperature (we use natural 
units for which h = ks = 1). For a given coupling 
strength r^, the decoherence rate is hence directly con- 
trolled by the temperature of the reservoir. The form of 
the equilibrium decoherence factor ([T]) has been success- 
fully verified in experiments involving mesoscopic super- 
positions of a radiation field inside a microwave cavity 
[5| and of motional states of a single trapped ion [6]. 
In most cases of practical interest, however, a quantum 
system is subjected to many noise sources with different 
amplitudes and correlation times, corresponding de facto 
to a nonequilibrium environment. 

In this paper we investigate nonequilibrium decoher- 
ence in the presence of many noise sources that act inde- 
pendently on the quantum system. The concrete sys- 
tem we have in mind is an ultracold ion confined in 



a large scale rf trap and subjected to different fluctu- 
ating electric fields [lo', |Tl[, but our approach is more 
general. We model each noise source by an equilibrium 
harmonic reservoir with given coupling constant, tem- 
perature and correlation time, and take non-Markovian 
effects fully into account. Despite the complexity of the 
problem, we show that the nonequilibrium decoherence 
factor can be expressed in terms of a simple measurable 
quantity. To accomplish this task, we use recent results 
from nonequilibrium statistical physics and map the cou- 
pling to M different equilibrium reservoirs to the coupling 
to a single nonequilibrium reservoir that we characterize 
with the help of a time-dependent effective temperature. 
We define the latter via a nonequilibrium extension of 
the fluctuation-dissipation relation and suggest an ex- 
perimental method to determine it in a linear Paul trap. 
The notion of effective temperature has been originally 
introduced in the description of glassy dynamics [8, 9]; it 
satisfies the expected properties of a temperature, such 
as a zeroth law (two interacting observables that evolve 
on the same time-scale have the same effective temper- 
ature) and has been measured in man y n onequilibrium 
systems, for instance in classical spin [12| and colloidal 
[13] glasses, as well as in granular systems Q and age- 
ing polymer glasses [l5|. We apply our results to the 
description of decoherence of a trapped ion interacting 
with two different amplitude reservoirs, treating the case 
of real and artificially engineered reservoirs [l6| . We show 
that in both instances nonequilibrium decoherence can 
be slowed down in a controlled manner as compared to 
the corresponding equilibrium situation, and explain how 
this behavior can be understood and predicted with the 
help of the effective temperature. 

The Hamiltonian of an oscillator S linearly coupled to 
M arbitrary harmonic reservoirs Rk is of the form 



M 



H = Hs 



EE 



2m 



kj 



Qkj 



(2) 

where Hs = /(2m) + mJ^x^ 12 is the Hamiltonian of 
the system oscillator with the usual notation and c/ej de- 
note coupling constants. Each equilibrium reservoir R]^ 
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is fully specified by a temperature and a spectral den- 
sity function Jk{oj) = (7r/2) c\- / {rrikji^kj) - ujkj) 
that characterizes the coupling to the system; we stress 
that we do not restrict ourselves to the Ohmic regime, 
J/e(u;) ~ cj, but include frequency-dependent damping. 
We describe the time evolution of the harmonic system 
with the exact quantum Langevin equation, 

mx{t)^ [ dt'r]{t-t')x{t')^muj^x{t) = F{t) -r]{t)x{0) , 



(3) 

which follows from the total Hamiltonian (j2j) [17]. Here 
the friction kernel r]{t) = ^]^rik{t) and the fluctuation 
force F{t) = ^]^Fk{t) are given by the sum of the re- 
spective contributions of the individual reservoirs. Under 
very general conditions (essentially that the reservoir is 
only weakly perturbed by the motion of the system [18]), 
any equilibrium reservoir can be modeled by an ensemble 
of harmonic oscillators [l9|. Since all reservoirs Rk are 
in thermal equilibrium, the correlation functions Ck{t) = 
{Fk{t)Fk{0) + Fk{0)Fk{t))/2 of the thermal forces and 
the friction kernels r]k{t) = 2 duo J k{oj) / {irmuj) co^ uot 
are connected via a fluctuation-dissipation relation. In 
the following, we consider the physically relevant regime 
for decoherence studies, in which temperature is larger 
than the oscillator frequency [H, [6|. The fluctuation- 
dissipation relation is then given in Laplace space by 
Ck[s\ = Tkrik[s\ [i^, indicating that, at equilibrium, en- 
ergy damping and thermal fluctuations evolve on the 
same time scale; we note that this equation can be in- 
terpreted as defining the equilibrium temperature of 
each reservoir. The above relation only holds at thermal 
equilibrium. However, away from equilibrium, we can 
extend the fluctuation-dissipation relation and define a 
frequency-dependent effective temperature as [2l|, 



T[s] = 



(4) 



where we have introduced the correlation function C[s\ = 
^j^Ck[s] of the total fluctuating force F[s]. We men- 
tion that a low-temperature extension of the effective 
temperature has been proposed in Ref. [22|. With the 



help of Eq. (|4]), we can now map the complex problem 
of a quantum system coupled to M arbitrary equilib- 
rium reservoirs to that of a system interacting with a 
single nonequilibrium reservoir characterized by a time- 
dependent effective temperature Teff(t) = Jq_ dt^T {t') 
(see Fig. d]). In the case of a single equilibrium reservoir, 
M = 1, we obtain T{t) = Tid{t) and the effective temper- 
ature reduces to the reservoir temperature, Teff = Ti, as 
expected. Conversely, we note that a general nonequilib- 
rium environment, with arbitrary r]{t) and C{t), can be 
approximated with M equilibrated reservoirs [21]. The 
effective temperature is a useful quantity to character- 
ize the violation of the fluctuation-dissipation relation in 
nonequilibrium systems. We will show below that it also 




FIG. 1: Schematic representation of the open quantum sys- 
tem. A complex environment, made of different equilibrium 
reservoirs, is mapped to a single nonequilibrium reservoir, de- 
scribed by an effective time-dependent temperature. 



provides deep insight into nonequilibrium decoherence. 

We turn to the dynamics of decoherence of a cat state 
in a general nonequilibrium environment. We focus on 
the regime of microscopic separation that corresponds 
to long coherence times, and investigate the influence of 
a nonequilibrium reservoir which evolves on similar or 
shorter time scales. The importance of nonequilibrium 
effects on short time scales can be made clear by consid- 
ering the coherent transport of ions in a segmented linear 
Paul trap: fast, nonadiabatic transport of the order of a 
few oscillation periods (i.e. a fraction of the coherence 
time) has been experimentally shown to drive ions into 
nonequilibrium excited states [23| . The accurate descrip- 
tion of short-time nonequilibrium effects appears there- 
fore crucial. Decoherence processes are conveniently ana- 
lyzed in Wigner phase space [3] . The Wigner function of 
a Gaussian cat state is given by the sum of two classical 
contributions W±{x^p) and a quantum interference term 

W,nt{x,p),W{x,p)=W- (X,p) + iy+(x,p)+iyint (X,p). 

During decoherence, the interference term progressively 
decays as a function of time. We quantify the disappear- 
ance of the contrast of the interference fringes by means 
of the peak-to-peak ratio between the interference and 
the classical terms of the Wigner function [24] . 



exp(-Aint) = 



peak 



2^yW- {x,p) |peakW+ {x,p) Ipeak 



(5) 



The nonequilibrium peak-to-peak ratio can be directly 
expressed as a function of the effective temperature (|4]) 
in the following way. We first note that for high tem- 
peratures, the dynamics of the system is dominated by 
momentum diffusion driven by the fluctuating force F{t). 
We can thus describe the evolution of the Wigner func- 
tion by the approximate Fokker-Planck equation, where 
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free evolution and dissipation terms are neglected 
dW{x,p,t) ... dW . ...d^W 



dt 



dW 
dp 



Here a^pit) = 2 dt' dt" C {f - t") is the momen- 
tum variance which describes the width of the Wigner 
function. The attenuation factor can now be read- 
ily evaluated from the solution of Eq. ([6]) and we find 
Ant(^) = d'^app{t)/2. According to the definition dU of 
the effective temperature, the noise correlator is further 
given by the convolution C{t) = dt'r]{t - t')T{t'). We 
therefore obtain, 

Ant(t) = d? f dt' I dt"r] it' - t") T {t") it - t') . (7) 



(6) 



The above nonequilibrium peak-to-peak ratio 
exp(— Aint) is the direct extension of the familiar 
equilibrium decoherence factor ([T]) to which it reduces 
when T{t) = T6{t); it completely characterizes de- 
coherence of a quantum system coupled to a general 
nonequilibrium environment, with arbitrary r]{t) and 
C{t). As in the equilibrium case, for a given coupling 
strength the decoherence rate is determined by the 
(effective) temperature of the nonequilibrium reservoir. 

Let us now illustrate our theoretical results by analyz- 
ing the concrete situation of a quantum oscillator that 
interacts with a fast and a slow reservoir. This is the sim- 
plest non-trivial out-of-equilibrium environment [21]. We 
begin by discussing 'real' environments before treating 
the case of 'engineered' environments in ion traps in the 
next section. We specify the fast (Markovian) reservoir 
by a temperature Tf and a delta-correlated friction ker- 
nel, r]f {t) = r]f5 {t)^ and the slow (non-Mar kovian) reser- 
voir by a temperature Tg and an exponentially-correlated 
kernel, r]s {t) = (r^s/r) exp (— t/r), with correlation time 
r. The effective temperature (|4]) is then given by 



T[s] 



Vf {sT + 1) + Vs 



(8) 



and the corresponding effective temperature in the time- 
domain reads, with r] = r]f -\- r]^^ 



Teff (t) 



Vs 



(T/-T,)exp 



T] t 

Vf r 



(9) 

The effective temperature (j9]) evolves monotonously from 
the fast reservoir temperature T/ at t = to the weighted 
average T^o = {rjfTf -\- rjsTs) /rj of the two tempera- 
tures at large times, t ^ {r]f/r])r. The exact peak-to- 
peak ratio exp(— 74int) (O for the nonequilibrium two- 
reservoir model is plotted in Fig. [2] for three differents 
sets of reservoir parameters. We observe that the pace 
of decoherence can be controlled by properly tuning the 
reservoir variables: nonequilibrium decoherence is slower 




FIG. 2: Interference- fringe contrast for a nonequilibrium two- 
reservoir system simulated by two random electric fields Ei (t) 
and E2{t). Nonequilibrium decoherence is slower (grey) or 
faster (black) than equilibrium decoherence (dashed) when 
the temperature Tf ^ {El) of the fast reservoir is higher or 
lower than the temperature Ts ^ (^|) of the slow reservoir. 
The dashed curve is the reference line corresponding to equal 
noise amplitudes {El) — {E2). The time-dependent effective 
temperature (shown in the inset) predicts this behavior: it 
is smaller than the reference temperature in the first case, 
while it is larger in the second. The effective temperature can 
be determined from measurements of the motional heating 
rate for trap frequencies u / {27t) varying between IMHz and 
lOOMHz. Parameters are: d — 7nm, a — 0.12nm, r = 5//s, 
(E!) = lO-^V^m-2 and (^|) = 6 • IQ-^^V^m-^ 



(faster) than equilibrium decoherence, which corresponds 
to Tf = T5, when the temperature of the slow reservoir 
is lower (higher) than the temperature of the fast reser- 
voir. We can therefore dramatically affect the decoher- 
ence properties of the system by, for instance, simply ex- 
changing the correlation times of the two reservoirs, while 
keeping their respective temperatures constant. This is 
a nontrivial non-Mar kovian effect. The effective tem- 
perature Teff offers an intuitive explanation for this be- 
havior. Following Eq. (j9j), we note that when Tf > Ts 
{Tf < Ts), the two-reservoir configuration is equivalent 
to a nonequilibrium reservoir with a decreasing (increas- 
ing) effective temperature. The approximate expression 
([7j) of the interference-fringe contrast then immediately 
shows that decoherence is slowed down (or accelerated). 

Linear Paul traps are promising candidates for quan- 
tum computation and an invaluable tool for the in- 
vestigation of the dynamics of decoherence [25]. Due 
to their weak interaction with the ambient reservoir, 
and the resulting long coherence time, they offer the 
unique possibility of exploring various decoherence sce- 
narios generated by the coupling to well-defined engi- 
neered reservoirs [l6|. The coupling to an equilibrium 
amplitude reservoir can thus be simulated with the help 
of a^uctuating electric field that acts on the charged 
ion [6]. In this case, the temperature of the reservoir 
is given by the variance of the field. The nonequilib- 
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rium two-reservoir configuration can be effectively sim- 
ulated with two random electric fields Ei{t) and E2{t) 
with corresponding correlation times. The correlation 
function of the superimposed noisy fields is given by 
C^{t) = e^{El)5(t) + (eVr)(E|)exp(-t/r) and is iden- 
tical to the two-reservoir correlation function C{t) with 
the identification e^{El) = r]fTf and e^(^2) = VsTg. 
Using the definition (jH of the effective temperature, 
we then find T^[s] = T + C^[s]/r]^[s] C^[s]/r]^[s], 
since the influence of the ambient reservoir can be ne- 
glected during the duration of the experiment. We note 
that the damping kernel r]^[s] = 77 is here constant, 
in contrast to the true two-reservoir case. As a result, 
the effective temperature for the two electric fields is 
Tfffit) = (eVr,) i{E^) + - exp(-t/T)), as shown 

in Fig. [2] (inset) for realistic trap parameters. Again, the 
decoherence properties can be understood with the help 
of the effective temperature: faster (slower) decoherence 
corresponds to higher (lower) effective temperature, as 
compared to the equal amplitude (Ef) = {E2) situation. 
The fact that T^{t) increases in the latter indicates that 
the noisy fields heat the system up, a direct consequence 
of the non-Markovian nature of the environment {T^ is 
a constant in the limit r ^ 0). We also mention that the 
exact ([5]) and approximate ([7j) peak-to-peak ratios are 
indistinguishable for the parameters of Fig. [2j We can 
now readily generalize these results to arbitrary complex 
situations: Given an unknown number of noise sources, 
with unknown amplitudes and correlation times, we can 
use the measurable effective temperature to predict the 
decoherence rate of the system via formula (O. 

An important property of the effective temperature is 
that it can be measured 12-15|. We next discuss an ex- 



perimental scheme that allows its determination in ion 
traps. We consider a single ion coupled to a nonequi- 
librium amplitude reservoir generated by an unspeci- 
fied number of electric fields with different noise ampli- 
tudes and correlation times (noisy electrical electrode po- 
tentials are the main sources of trap decoherence [11]). 
The corresponding effective temperature ^ is given by 
the ratio of the noise correlation function and the fric- 
tion kernel. Both quantities can be determined experi- 
mentally by measuring the motional heating rate of the 
trap [26]. The latter is related to the power spectrum 
of the fluctuation force, S{uj) = J dt ex.-p{iujt)C{t)^ by 
h S {uo) / {Amhw) ^ where e is the charge of the ion 
[2^ . The noise correlation function can hence be directly 
obtained from measurements of the heating rate for dif- 
ferent trap frequencies uo. A first set of measurements 
of the heating rate for the ambient equilibrium reservoir 
alone yields the friction kernel 77 (which in this case has 
been shown to be frequency independent |26|]) via the 
fluctuation-dissipation relation, S{(jj) = 2D = 2r]T. A 
second set of measurements of the heating rate in the 
presence of the random electric fields further provides 
the correlation function C{t) of the nonequilibrium engi- 



neered reservoir and, in turn, the time-dependent effec- 
tive temperature. 

To conclude, the present findings show that nonequilib- 
rium decoherence induced by an unknown number of dif- 
ferent noise sources can be successfully described within 
a unified framework, in a way similar to equilibrium deco- 
herence. They further highlight the central role of the ef- 
fective temperature to characterize, understand and pre- 
dict the loss of coherence of a quantum system in a gen- 
eral nonequilibrium environment. 
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